Quiescence: a mechanism for escaping the effects of drug on cell 

populations 

Tomas AlarcorJ*) 
Basque Centre for Applied Mathematics, 
Bizkaia Technology Park, 48160 Derio, Spain 

Henrik Jeldtoft JenserJl) 

Institute for Mathematical Sciences, Imperial College London, 

53 Princes Gate, London SW7 2PG, 
UK & Department of Mathematics, Imperial College London, 
South Kensington Campus, London Sw7 2AZ, UK 
(Dated: February 24, 2010) 

Abstract 

We point out that a simple and generic strategy to lower the risk for extinction consists in 
the developing a dormant stage in which the organism is unable to multiply but may die. The 
dormant organism is protected against the poisonous environment. The result is to increase the 
survival probability of the entire population by introducing a type of zero reproductive fitness. 
This is possible, because the reservoir of dormant individuals act as a buffer that can cushion fatal 
fluctuations in the number of births and deaths which without the dormant population would have 
driven the entire population to extinction. 
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FIG. 1: Schematic representation of the population dynamics considered in this paper. The arrows 
between the different types of individuals represent the possible population flows between them. 
The arrows pointing to the types themselves represent proliferation. The other arrows stand for 
death of each type of individuals. The thickness of the arrows stand for the realtive magnitude of 
each of the flows or processes. 

Introduction. Drug resistance is a very serious problem not least for chemotheaputic 
treatment of cancer and it is very important to understand, as far as possible, ways to 
countermeasure drug resistance in a number of biomedical contexts . In order to unravel 
the possible mechanisms responsible for this phenomena Iwasa and coworkers developed a 
prominent theoretical model of the evolutionary dynamics of escape [l, 2|. They base their 
approach on the assumption that n point mutations in some crucial parts of the genome 
are necessary for escape. They further assumed that the different mutants can be described 
by binary strings (with entries +1 or —1) of length n. There are 2 n — 1 such mutants. 
It is assumed that treatment reduces the proliferation ratios of sensitive mutants: R < 1 
whereas resistant mutants are such that R > 1. The corresponding evolutionary dynamics 
is modelled in terms of Galton- Watson multitype branching process (GWMBP) |3|] where 
at each generation each individual of each type has a given (in general, mutant-dependent) 
probability of mutating and producing offspring belonging to a different type. The problem 
is to calculate the probability that a resistant mutant is reached within a population of size 
N. The model proposed by Iwasa et al. has been analysed in more detail in 0, 5|. Whilst 
this approach is both biologically and mathematically sound, we believe the efficiency with 
which evolutionary escape allows organisms to avoid extinction, when attacked by lethal 
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drug, can be understood in a more simple way as a consequence of a dormant phase. To 
assess the generic importance of dormancy as an escape strategy we develop a very simple 
model framework, in which we only focus on the most essential aspects of population and 
evolution dynamics. 

Here we show that the existence of a dormant phase can be crucial for a population to 
escape extinction. We consider the following simplified scenario. We study the survival 
probability of a population of organisms that can exist in the form of three different types, 
see the diagram [IJ The types differ in their response to the presence of a drug. Type (1) 
and type (2) have similar reproduction and death rates when no drug is present. However, 
the drug is supposed to be lethal to type (2) but neutral to type (1). We are interested in 
how the existence of a dormant mode, type (3), effects the survival probability of the entire 
population. The dormant type cannot reproduce, nor is it susceptible to the drug. However 
type (3) can die and it can undergo a transformation back to type (1). To understand the 
effect of the dormant type we consider different realisations (indicated in the diagram [T|) of 
the possible flow between the three different types. 

Our model is relevant to a number of biological cases, in particular to entities such as 
cancer cells, which have been observed to evolve resistance to therapy: Treatments impose 
a selective pressure which eliminates the lesser fit strands of the corresponding populations 
but also drives an evolutionary process whereby better adapted individuals, i.e. individuals 
immune to the effects of the corresponding drug, eventually take over. Further rationale for 
our model is provided by the response of tumour cells to hypoxia (i.e. oxygen starvation). 
It is a well-known fact that hypoxia induces arrest of the cell-cycle js, 7]. This means that 
the rate at which cells replicate under hypoxia is drastically reduced, as hypoxia down- 
regulates the activity of the pathway regulating the progression through the cell-cycle {(J. 
Moreover, another well-understood fact about cell survival/death regulation is that there 
exists cross-talk between the pathways regulating cell death and cell division @-0| • This 
means that, in normal circumstances, the down- regulation of the cell cycle machinery implies 
the down-regulation of the apoptotic (programmed cell death) machinery. It can be argued 
that such cross-talk is disrupted in cancer cells, and that cancer cells are such that cell death 
regulatory pathways are down-regulated. Therefore, under hypoxic conditions, cancer cells 
undergo a drastic reduction of their division and death rates. 

The study of the influence of hypoxic cancer cells fits rather well within the original remit 
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in which the dynamics of evolutionary escape was put forward, as hypoxic cells are known 



to play a major role in the resistance to chemo- and radio-therapy in tumours [111]. To study 
the influence of such a sub-population on the global dynamics of the total population, we 
model it as a dormant or quiescent population with neglegible proliferation rate and small 
death rate. 

Previous works have analysed some of the effects of quiescent cells in tumour growth (see 



for example 



12 



13]). All these previous studies hint to the role of quiescence cells in the 



dynamics of tumour growth and its role in resistance to tumour growth but the issue of the 
evolutionary dynamics involved is not directly addressed. 

We now turn to the study of the survival probability. We disentangle the interplay 
between drug susceptibility, type (2), and dormancy, type (3), by analysing three different 
versions of the population dynamics, all are depicted in the diagram [1] and formulated as 
the following three Models. 



Model A — This is the "normal" situation, where no drug is present. Type (1) and type 
(2) are essentially equivalent, except that when type (1) reproduces it may undergo a 
"mutation" and end up as type (3). When type (2) reproduce it may mutate and become 
type (1). 



Model B — This is the situation in the presence of the drug. The death rate of type (2) 
is now significantly bigger than the death rates of type (1) and type (3). Moreover, the 
drug makes type (2) unable to reproduce and therefore the flow from type (2) to type (1) is 
absent. 

We will below find that the dormant stage increases the survival probability of the 
population in the presence of the drug. To emphasis that the increased ability of the 
population to escape extinction is in fact caused by the presence of a dormant stage, we also 
consider an extreme version of the dynamics in which type (1) is unable to die. Namely, 



Model C — This version of the the dynamics is equivalent to model B except that type (1) 
now is assumed not to be able to die directly, but has to flow through either type (2) or type 
(3) to do so. We then demonstrate that even in this extreme situation does the availability 
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of the dormant stage enhance the population's chance for avoiding extinction. 

We notice that the different types (1), (2) and (3) can be thought in epidemiological terms 
as susceptibles, infected and immune. The version of the dynamics defined as model B can be 
thought of as representing an age structured population. Type (3) is then juveniles, type (1) 
mature reproduction active individuals and type (2) are individuals in the post-reproductive 
stage. 

An economical and precise way to present the dynamics is in terms of the generator 
functions for the corresponding Galton- Watson multitype branching process, we include 
these generator functions in the tables in the appendix. Within the theory of multi-type 
branching process, the condition for eventual non-extinction is given in terms of the spectral 
radius, p, of the matrix A = (a^), whose entries are the expected values of the number of 
offspring of type j produced of an indvidual of type i. If p > 1, there is a finite probability 
of N(t — » oo) > 0. The quantities are calculate in terms of the corresponding generating 
functions (Table HT1): Aij = djGi\g=i. From these standard arguments of the theory of multi- 
type branching process, which carry on in a straightforward manner when size-dependence 
is taken into account we can establish from Tables HI and HTl Appendix A, that the 
condition for asymptotic survival, that is for P[N(t — > oo) > 0] > 0, is 



2(1 -r 12 )e-» N °° > 1, (1) 

i.e. ri2 < 1/2. Otherwise, the probability of eventual extinction is 1. Consider now the 
presence of a third type (type 3) within the population, namely, quiescent individuals. 
These individuals are not allowed to proliferate but are resilient to the environment and 
can survive under hostile conditions, hence we assume that their death rate e << 1. In 
addition, quiescent cells are assumed to be able to revert back to type 1 at a given rate 
which depends on the availabity of resources. The issue we intend to analyse is whether 
the introduction of a quiescent sub-population helps to escape from the whole population 
being extinct. More precisely, the question we aim to address is: Assume r\2 > 1/2, is it 
possible for a population whose dynamics is described in the by Model B and Table [III to 
elude eventual extinction? Specifically this scenario has been considered within the context 
of modelling of tumour growth, in particular in the response of cancer cells to hypoxia (low 
levels of oygen) [12J. Cancer cells appear to become quiescent in response to hypoxia, which 
is thought to give them an advatage in their competition with their normal counterparts as 
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FIG. 2: Simulation results for r 12 = 0.5. Black circles correspond to r 13 = 0.5 and red squares to 
r 13 = 0.25. Other parameter values: /i = 0.001, e = 10 -5 , and a = 1. 

well as resistance to radio- and chemo-therpay [ijj]. 

We now present the results of simulations of the survival probability for the three different 
scenarios represented by model A, B and C. 

Simulation results. Numerical simulations of the underlying multi-type branching process 
(see Appendix A for details) confirm that, in fact, the presence of a quiescent population as 
a mechanism for escape to harsh conditions is indeed feasible. Fig. [2] shows how the survival 
probability depends on the probability of an indvidual of type (1) to become quiescent, 
ri3. We can see that as ri 3 decreases the threshold for survival in terms of r 3 i, i.e. the 
probability of a quiescent individual to revert to type (1), moves towards smaller values of 
the parameter r 31 . This means that, following a decrease in the flux of individuals from 
the type (1) population to type (3), survival is only possible by reducing the inverse flux. 
This observation reveals that the mechanism by which quiescence helps escape is by acting 
as a reservoir where part of the population can be safely "stored". If the flux from type 
(3) back to type (1) is too big, it effectively increases the flux from type (1) to type (2), 
thus increasing lethality. This behaviour, however, is sensitive to the decay rate of type (3) 
individuals, e: as this parameter increases the survival probability decreases (see Fig. H|). 
This means that the death rate of quiescent individuals being small is an instrumental factor 
for quiescence-induced escape from harsh environments. 

Fig. |3] shows that as ryi increases further into the regime where Eq. ([T]) predicts sure 
extinction of the (2)-type population, survival is only possible by decreasing r 3 i, i.e. by 
increasing the average time individuals spent in the quiescent state. 

To test further this scenario, we perform simulations in which the following situation is 
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FIG. 3: Simulation results for ri3 = 0.5. Black circles correspond to r\i = 0.5, red squares to 
ryi = 0.525, and green diamonds to r\2 = 0.55. Other parameter values: fi = 0.001, e = 10 -5 , and 
a = 1. 




FIG. 4: Simulation results for = 0.5 and r\2 = 0.5, and green diamonds to r\2 = 0.5. Black 
circles correspond to r-^i = 0.05 and red squares r$i = 0.1. Other parameter values: /i = 0.001, 
and a = 1. 

considered. First we let a population whose dynamics is given by Model A, evolve for some 
time in an environment free of the hostile agent. In this environment, both type (1) and 
type (2) individuals can thrive. The dynamics in an environment free of any hostile agent 
is described in terms of the generating functions shown in Table [B The introduction of the 
hostile agent is describe by changing the dynamics of the population to the one described by 
model B (see the diagram [1] and Table HT|) . which is essentially the same as model A (Table HJ 
but with type (2) individuals doomed to perish. We further assume that the hostile agent is 
active only for a given period of time, after which its detrimental effect on type (2) organisms 
ceases and the dynamics of the population reverts to Model A (as per Tabled]). Parameter 
values are chosen so the dynamics of the agent-free population is super-critical (i.e. the 
corresponding survival probability is bigger than zero). After the hostile agent has been 
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FIG. 5: Simulation results showing one realisation of the population dynamics of the reaction to the 
(transient) application of a therapeutic agent. We can see that in plot (a), which corresponds to a 
population with no quiescent state, the population upon application of the drug. On the contrary, 
a population with quiescence is able to survive to the period over which the drug is active, and 
recovering once the effect of the drug has worn off (plot (b)). Black lines correspond to the total 
population, red lines to individuals of type (1), green lines to individuals of type (2), and blue lines 
to quiescent (type (3)) individuals. Parameter values: T\% = = 0.5, r\3 = = for the upper 
pannel and 7*13 = r%i = 0.25 for the lower pannel, r^i = 0.05, // = 0.001, e = 10~ 5 , and a = 1. 
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FIG. 6: Simulation results showing the survival probability, P[N > 0], for the population dynamics 
of the reaction to the (transient) application of a therapeutic agent as a function of the parameter 
ri3, i.e. the parameter controlling the flux of population between types 1 and 3. Parameter values: 
r 12 = r 2i = 0.5, r 3 i = 0.05, n = 0.001, e = 10" 5 , and a = 1. 
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FIG. 7: Simulation results for model C with r\i = 0.51. Black circles correspond to r\3 = 0.45 and 
red squares to r%s = 0.25. Other parameter values: /i = 0.001, e = 10~ 5 , and a = 1. 

removed we calculate the probability that the population survives both with and without 
quiescence. This scenario is highly relevant to evolutionary escape from drugs which have a 
finite life-span. This particular escape problem was considered Iwasa et al.(l|. 

The corresponding results are shown in Figs. |5]and[6j Fig. [6] shows the survival probabil- 
ity, P[N > 0], as a function of r 13 and, indeed, we observe that when quiescent individuals 
are not considered, i.e. r 13 = 0, the probability of surviving after the end of the activity 
period of the hostile agent is null. As r\s increases, so does the corresponding survival prob- 
ability, thus further confirming that quiescence is a feasible mechanism to help biological 
populations to escape harsh environmental conditions. Moreover, Fig. which show par- 
ticular realisations of the process with (panel (b)), and without (panel (a)), quiescence also 
helps to understand more about the mechanism whereby quiescence allows escape: as can 
be seen in Fig. [5j panel (b), during the period of activity of the hostile agent most of the 
population is actually quiescent with a little proportion of the population in types (1) and 
(2). This means that quiescence mediates survival by acting as a reservoir or buffer. 

An alternative model. Whilst the results discussed until here provide compelling evidence 
in favour of hypothesis that the presence of quiescence can induce escape without the need 
of an increase in the reproductive fitness of any particular type of organism, some doubt 
could be cast on our argument so far. One might feel inclined to argue that our model 
still relays on a multi-type population where one type (type 1) is effectively more fit than 
the others. To adress this possibility, we consider model C in diagram [TJ This version is 
closely related to the current discussion concerning the response of cancer cells to oxygen 
starvation. This model is also formulated in terms of a GWMBP and characterised in terms 
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of the generating functions given in Table IHIl Its rationale is as follows. Let us assume that 
a population of cells is divided in two types: those in which the pathways regulating cell 
division are activated (type 1) and those in which the active pathways are those regulating 
apoptosis, i.e. cell death, (type 2). Type 1 cells can either proliferate, or stay as they are, 
or suffering activation of the apoptotic pathways, thus becoming type (2) cells. The flux of 
population between types (1) and (2) is controlled by the parameter r±2 ■ Type 2 cells are 
those marked for cell death, but the model provides for their staying within the population 
for a while (not dying immediately), but they cannot proliferate. Type (3) cells are, as in 
the previous model, quiescent cells. In this context the introduction of a hostile agent (drug, 
removal of oxygen, etc.) corresponds simply to increasing the value of r 12 . 

In the absence of quiescence, standard arguments reveal that extinction with probability 
1 occurs when ru > 1/2. The question is once again whether quiescence can rescue the 
population from extinction under such conditions. Fig. [7] shows that this is indeed the case, 
provided that the flux out of the quiescent state back into type (1) does not exceed a critical 
value. This results are completely analogous to those obtained for the previous model. 

Invasion dynamics. We know discuss under which conditions a small proportion of in- 
dividuals, y, which can undergo quiescence, modelled by Model A and B depending on 
whether a drug is present or not, respectively, and with ^ 0, can take over a population 
of fully growing individuals, i.e. a population modelled by Model A and B with r 13 = 0. In 
particular, we will analyse the corresponding invasion probability, Pp{y), by direct simula- 
tion of the population dynamics and by using the analytical approximation provided by the 
so-called evolutionary formalism 16Ml9| which allows a more thorough exploration of the 
behaviour of the invasion probability as a function of the model parameters. 

The problem we address here relates to whether a small population of mutants with can 
take over an incumbent population. Demetrius and coworkers 17|, ll8j have developed a 
brmalism based on the application of ideas from ergodic theory to evolutionary problems 
Ji| and the diffusion approximation [^(J. This formalism allows us to estimate the fixation 
probability of a mutant population in the presence of an incumbent species. Here we gener- 
alise this formalism to apply it to the problem of whether a genetic inactivation generating 
new phenotypes can invade the incumbent population. 

The starting point of this formalism is the following fundamental equation: 
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r = H + F, (2) 

where r = log(Ao) is the growth rate (or Malthusian parameter), with Ao is the dominant 
eigenvalue of A, whereas H and F are the entropy and the proliferative potential, defined 
as: 



where r>,e is defined as: 



H = -^7r^lnp rf 



i.t 



1,1 



(3) 



Pv 



(4) 



with VA = AqV, and rr is the stationary distribution associated to V — (pij): 7T ? ; — ViUi, with 
U given by AU = AqU. Eqs. (|2])-(j4|) are derived from a variational principle 16], namely 



sup ( H u 



(5) 



where v is a Markov measure H v — — YTj=i KA?) ^sKA/) aric ^ 7 = ^°S a x xi, where Aj, j = 
1, . . . , n is a partition of the phase space of the system and a XQXl is the transition probability 
between two states of the system, xq and X\. The solution to this variational problem 



produces Eqs. ()2])-( 
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According to Demetrius et al. 18] the diffusion approximation yields the following equa- 
tion for the fixation or invasion probability as a function of the initial concentration of 
mutants, y is given by: 



Pf(v) 



1-1 



2W- 

' M 



+1 



Act 2 



+ 1 



(6) 



^l 1 

where the total population A" is assumed to be constant and s is defined by: 



Ar - 



Aa 2 

(AO 



(7) 
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where Ar = r q — r, Aa 2 = a 2 — a 2 , with r q (r) is the growth rate of the quiescence (normal) 
population, and a 2 (a 2 ) is the variance of the quiescent (normal) population, (N) is the 
average stationary population. 

The demographic parameters (i.e. r, a 2 , r q , and a 2 ) can be estimated from the evolu- 
tionary formalism (see 18|, ll9j). The growth rates are given by r = lo g An , where An is the 
dominant eigenvalue of the of the matrix A = (cijj) = (djGi(x)\$=i). In [18J, it is shown that 
the parameter a 2 can be obtained by slightly perturbing the parameters that determine the 
dynamics of the system, i.e. the mean-field dynamics being given by A(S) = (a]? 6 ), and 



then doing an expansion for small 5. Accordingly, a 2 is given by: 

2 _ dH(5) 



a' 



(8) 

5=0 



dS 

Hence, in the linear approximation with H(5) given by H(5) = H + 8H$, we have that 
cr 2 (7) = —Hs, which is given by: 

H& = - W 1 + lo S Pij) + Kipfj log Pij (9) 

where 7Tj(5) = iXi +87rf ) and Pij(S) = +Sp\j are the corresponding linear approximations 
to 7r and P = (p^) when we take A(5) = 5 « 1. The details of how these quantities 

are actually calculated is given in Appendix A of |19| . 

In order to show that this formalism can be used to analyse the invasion of a population 
possessing a dormant type, we have done simulations where one population, which is sensitive 
to the drug, competes with two different population: one capable of undergoing quiescence 
and another one which is unable to enter such state. We compare the numerics with the 
analytical results obtained from Eq. (jSj) in Fig. [HJ Both analytical and numerical results 
show that whereas in the former case the quiescent population takes over the sensitive 
population almost surely, in the latter case is unlikely that the non-resident population takes 
over the resident one. It is worth remarking that the good agreement between simulations 
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[8] implies that the evolutionary formalism developed by 
18] is able to make the invasion problem considered here 



and analytical result seen in 
Demetrius and coworkers 
analytically tractable. 

Summary and Discussion. We have shown that quiescence is a feasible mechanism for 
biological populations to escape hostile environments. This mechanism is expected to be 
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FIG. 8: Fixation probability for the competition between two populations in the presence of drug. 
Squares correspond to the competition with a population capable of undergoing quiescence whereas 
circles correspond to the competition with a population that cannot undergo quiescence/ Solid line 
corresponds to the analytical expression of Pf{v) for the latter case 




FIG. 9: This plot illustrates the phenomenon of persistence. After introducing a drug (typically, 
an antibiotic), there is an initial period of fast decay of the population followed by a cross-over 
after which the rate of population decays is much slower. In bacterial populations such as E. Coli, 
this is thought to be due to a phenotype switch: Cells switch to a slow growing phenotype that 
is more resilient to the effect of the grug than the normal phenotype. This phenomenon cannot 
be of genetic origin, as cells regrown from the persistent population upon removal of the drug are 



sensitive to it 



22|. 



very relevant to the important issue concerning population dynamics of cancer cells and the 
competition mechanisms between cancer cells and their normal counterparts. In particular 
our model addresses issues related to the response of cancer cells to oxygen starvation 

HQ. 

The mechanism involved in quiescence-dependent escape is in essence simple: it provides 
a buffer for the population to be safe from the hostile environment. I.e. in cases where 
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the population would go extinct due to a fluctuation in the birth and death events the 
reservoir in the quiescence buffer population can bring the population back from the brink 
of extinction. The efficiency of this mechanism for escape is sensitive to the entry and exit 
rates to and from the quiescent state, respectively. 

The main difference between our mechanism for escape and the one proposed in 

HQ 

is the following. Iwasa and co-workers propose a mechanism based on random search of 
the state of types of individuals for one that is fitter than the others in the presence of a 
given selection pressure (drug, etc.). In contrast we have demonstrated here that the overall 
survival probability of a population may increase by introducing a dormant type. Since this 
type is unable to reproduce its reproductive fitness is zero, nevertheless the existence of this 
drug resistant stage is able to improve the fitness of the entire population in as much as the 
population obtains a higher probability for survival. 

The mechanism proposed here shares a number of features in common with the phe- 



nomenon of bacterial persistence 22l-l25|. Persistence is a form of resistance to antibiotics 



exhibited by bacterial colonies where resistance is not acquired by a gene mutation which 
allows the bacteria to grow exponentially fast in the presence of a particular drug. Instead, 
as schematically shown in Fig. |9l killing of bacteria goes through a fast phase where the 
population decreases exponentially until the killing rate slows down leaving behind a rem- 
nant of cells which differ from the sensitive ones in that they are in a dormant state but 
are otherwise genetically identical to their sensitive counterparts. In fact, when the drug is 
removed, a colony of "normal" bacteria is regrown from the persistent cells. The mechanism 
for persistence appears to involve a phenotypic switch 22j where cells switch from a rapidly 
growing, but drug-sensitive phenotype to a dormant but drug-resistant phenotype, although 
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there has recently been the suggestion that persistence might be a social trait 

According to 22| , there exist two types of persisters. Type I persisters, according to their 
terminology, are only produced during the exponential growth phase and therefore their 
numbers are fixed at the time of innoculation and determined by the size of the innocolum. 
Type II persisters, on the contrary, divide and grow continuously, but an order of magnitude 
slower than their non-persistent counterparts, and their numbers are determined by the total 
population numbers. The quiescence mechanism put forward here is therefore distinct from 
type I persisters, but very similar to the behaviour exhibited by type II persisters, which 
means that the analysis methods, in particular, the evolutionary formalism used to study 
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the competition between sub-populations exhibiting different behaviour, can be extended to 
study this type of bacterial persistence. 

Another issue in which our model differs from previous work on bacterial persistence is 
the following. We are consider different populations adopting different strategies, namely, a 
wild-type population which thrives in the absence of drug composed of two subpopulations: 
one that is capable of undergoing quiescence and another one that is not. By adopting this 
scenario we can study the evolutionary dynamics of these populations and analyse which 

n that respect, 
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strategy is evolutionary stable and which one is susceptibel to be invaded, 
our analysis goes beyond, for example, the population models proposed by 

Within the context of the problem of how the hypoxic sub-population affects the dynamics 
of the whole tumour, our model sheds further light on this issue. It is commonly thought 
that hypoxia increases the probability of survival of the tumour by providing a selective 
pressure that favours the evolution and survival of more aggresive phenotypes ll|. Our 
model shows that quiescence by itself is enough to increase the survival probability of the 
population without further increase in the phenotypic variety of the population. Obviously, 
the residual population provides a springboard for these evolutionary processes to ensue. 

TA and HJJ gratefully acknowledge the EPSRC for funding under grant EP/D051223. 



Appendix A: Simulations 

In this appendix we briefly describe the method we have used to produce our simulation 
results. Our simulations start with one single individual of type 1. Its offspring is the 
determined by the probabilities of producing descendants as prescribed by the generating 
functions given in Tables H| \U\ and IHH corresponding to Models A, B, and C, respectively. 
In general, the coefficients of the Taylor expansion of Gi{x) = J2i n m=o Pi;inmx[x2 x T are the 
probabilities per generation per indvidual (of type i) that an type i-individual produces / 
descendants of type 1, n descendants of type 2, and m descendants of type 3. In subsequent 
generations, we go over all the individuals in existance in the last generation and the numbers 
and types of their descendants within the next generation are calculated in the same way. 
Simulations are run over 1000 realisations of 1000 generations each. The survival probability 
P[N > 0] is calculated as the ratio between the number of simulations such that N(T = 
1000) > and the total number of realisations. 
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TABLE I: Model A. Generating functions of the probabilities of the number of offspring for 
individuals belonging to our population previous to the introduction of a hostile agent. In this 
situation, both type 1 and type 2 indviduals are well adapated and can proliferate subject to 
the restrictions imposed by the carrying capacity. Type 3 correspond to quiescent individuals. 
The new parameters introduced here are r 2 i, i.e. the mutation probability of type 2 into type 1, 
the probability of a type 2 individual to enter quiescence, r 23 , and the probability of a quiescent 
individual to revert to type 2. The rest of the parameters have the same physical meaning as in 
Table O 
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TABLE II: Model B. Generating functions of the probabilities of the number of offspring for 
individuals of type 1, i.e. well adapated to harsh conditions, Gi(x), type 2, i.e. those for which 
the environment is lethal, G 2 (x), and type 3, i.e. quiescent individuals, G 3 (x). The parameter 
fi is the carrying capacity, which accounts for the limitation in resources, a is a measure of the 
remanent resilience of inviduals of type 2 to the environment, and e is the death rate (probability 
of death per individual per generation) of the quiescent cells. The quantity ru with j = 1,2 can 
be interpreted as the muation rate of individuals of type 1 into individuals of type j. The quantity 
r 3 i is the probability per quiescent individual per generation of a type 3 individual to revert to ype 



1. 


Population 1 


G x {x) = 1 - e~^ N + e-» N r X 2xl + e~^ N r l3 xl + e~» N {l - r 12 - r 13 )x\ 


Population 2 


G 2 (x) = a + (1 - a)x 2 


Population 3 


G 3 (x) = e + r 3 ie'^ N xi + (1 - e - r 3 ie~ flN )x 3 
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TABLE III: Generating functions corresponding to model C (see main text for details). Type 1 and 
type 2 correspond to the proliferative and apoptotic phenotypes, respectively. Type 3 corresponds 
to the quiescence population. 
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